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- - - Treating the production of electron and positron pairs by a strong electric field from the 
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^^ ' vacuum as a quantum tunneling process we derive, in semiclassical approximation, a general 
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^\J , expression for the pair production rate in a z-dependent electric field E{z) pointing in the 
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^-direction. We also allow for a smoothly varying magnetic field parallel to E(z). The result 
is applied to a confined field E{z) ^ for \z\ ^ £, a, semi-confined field E{z) ^ for z ^ 0, 
and a linearly increasing field E{z) ~ z. The boundary effects of the confined fields on 
pair-production rates are exhibited. A simple variable change in all formulas leads to results 
for electric fields depending on time rather than space. 

In addition, we discuss tunneling processes in which empty atomic bound states are spon- 
taneously filled by negative-energy electrons from the vacuum under positron emission. In 
particular, we calculate the rate at which the atomic levels of a bare nucleus of finite size r^ 



>, 
0^ , and large Z 3> 1 are filled by spontaneous pair creation 
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INTRODUCTION 



The creation of electron-positron pairs from the vacuum by an external uniform electric field in 
j^ ■ spacetime was first studied by Sauter jilj as a quantum tunneling process. Heisenberg and Euler [4] 

extended his result by calculating an effective Lagrangian from the Dirac theory for electrons in 



a constant electromagnetic field. A more elegant quantum field theoretic reformulation was given 
by Schwinger [3|] based on calculations of who calculated within Quantum Electrodynamics (QED) 
the one-loop effective action in a constant electromagnetic field. A detailed review and relevant 



a 



references can be found in Refs. [J] and 

Apart from its purely theoretic interest, the pair-production in nonuniform fields is experimen- 
tally relevant in collisions of laser beams [a] and heavy ions [71,13], and as a possible explanation of 
the powerful Gamma Ray Bursts in astrophysics [a, l9[]. It is also important for understanding the 
plasma oscillations [lO| of electrons and positrons after their creation in electric fields. 

The rate of pair-production may be split into an exponential and a pre-exponential factor. The 



exponent is determined by the classical trajectory of the tunneling particle in imaginary time which 
has the smallest action. It plays the same role as the activation energy in a Boltzmann factor with 
a "temperature" h. The pre-exponential factor is determined by the quantum fluctuations of the 
path around that trajectory. At the semiclassical level, the latter is obtained from the functional 
determinant of the quadratic fluctuations. It can be calculated in closed form only for a few classical 
paths [Ul]. An efficient technique for doing this is based on the WKB wave functions, another on 
solving the Heisenberg equations of motion for the position operator in the external field [11[. 

Given the difficulties in calculating the pre-exponential factor, only a few nonuniform electric 
fields in space or in time have led to analytic results for the pair-production rate: (i) the electric field 
in the z-direction is confined in the space x < xq, i.e., E = E{x)z where E{x) = Eq9{xq — x) [I4]; 
(ii) the electric field in the z-direction depends only on the light-cone coordinate z+ = (t + z)/\/2, 
i.e., E = E{z^)z [13]. If the nonuniform field has the form E[z) = ii^o/cosh (z), which we shall 
refer to as a Sauter field, the rate was calculated by solving the Dirac equation [14] in the same way 
as Heisenberg and Euler did for the constant electric field. For general space and time dependences, 
only the exponential factor can be written down easily — the fiuctuation factor is usually hard to 
calculate |2|. In the Coulomb field of heavy nucleus whose size is finite and charge Z is supercritical, 
the problem becomes even more difficult for bound states being involved in pair production, and a 
lot of effort has been spent on this issue 7|, la. llSl]. 

If the electric field has only a time dependence E = E{t), both exponential and pre-exponential 
factors were approximately computed by Brezin and Itzykson using WKB methods for the purely 



periodic field E{t) = E^qcoswo* [16]- The result was generalized by Popov in Ref. [17| to more 
general time-dependent fields E{t). After this, several time-independent but space-dependent fields 
were treated, for instance an electric field between two conducting plates [118], and an electric field 
around a Reissner-Nordstrom black hole [19[]. 

The semicalssical expansion was carried beyond the WKB approximation by calculating higher- 
order corrections in powers of h in Refs. [201] and 21]. Unfortunately, these terms do not comprise 
all corrections of the same orders h as will be explained at the end of Section IIIBi 

An alternative approach to the same problems was recently proposed by using the worldline 



formalism 22] , sometimes called the "string-inspired formalism" . This formalism is closely related 
to Schwinger's quantum field theoretic treatment of the tunneling problem, where the evaluation 
of a fluctuation determinant is required involving the fields of the particle pairs created from 
the vacuum. The worldline approach is special technique for calculating precisely this functional 
determinant. Within the worldline formalism, Dunne and Schubert 23|] calculated the exponential 



factor and Dunne et al. [2^ gave the associated prefactor for various field configurations: for 
instance a spatially uniform, and single-pulse field with a temporal Sauter shape oc 1/ cosh cot. 
z-dependences, a numerical calculation scheme was proposed in Ref. [251] and applied 
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further in 2^ . For a multidimensional extension of the techniques see Ref. |27l ] . 

In this article we derive a general expression for the pair-production rate in nonuniform electric 
fields E(z) pointing in the z-direction and varying only along this direction. A simple variable 
change in all formulas leads to results for electric fields depending on time rather than space. As 
examples, we shall treat three cases: (i) a nonzero electric field confined to a region of size i, 
i.e., E{z) 7^ 0, \z\ ^ i (Sauter field see Eq. (J56l) ): (ii) a nonzero electric field in a half-space, i.e., 
E{z) ^ 0,z ^ (see Eq. (j77p ): (iii) an electric field increasing linearly like E{z) '^ z. In addition we 
shall study the process of negative-energy electrons tunneling into the bound states of an electric 
potential with the emission of positrons. We consider two cases: (1) the electric field E{z) ~ z of 
harmonic potential V{z) ~ z^; (2) the radial Coulomb field E{r) = eZ/r'^ with large Z outside the 
nuclear radius r^. 

II. SEMI-CLASSICAL DESCRIPTION OF PAIR PRODUCTION 



The phenomenon of pair production in an external 
mechanical tunneling process of Dirac electrons 23, 



electric field can be understood as a quantum- 



29|. In the original Dirac picture, the electric 



field bends the positive and negative-energy levels of the Hamiltonian, leading to a level-crossing 
and a tunneling of the electrons in the negative-energy band to the positive-energy band. Let the 
field vector 'E(z) point in the z-direction. In the one-dimensional potential energy 



V{z) = -eAo{z) = e dz'E{z') (1) 

of an electron of charge — e, the classical positive and negative-energy spectra are 



£±{Pz,P±; z) = ±^{cp-,f + c^p\ + (mec2)2 + V{z), (2) 

where Pz is the momentum in the z-direction, p_L the momentum orthogonal to it, and p±_ = |p_l|. 
For a given energy £, the tunneling takes place from z^ to 2;+ determined hy pz = ^ in Eq. ([2]) 

S = S+{Q,pr,z+) = S^{Q,pr,z^). (3) 

The points z± are the turning points of the classical trajectories crossing from the positive-energy 
band to the negative one at energy £. They satisfy the equations 



V{z±) = ^^c^p\ + mlc^ + E. (4) 



This energy-level crossing £ is shown in Fig. [T] for the Sauter potential V{z) oc iscnh{z/t). 




FIG. 1: Positive- and negative-energy spectra £±{z) of Eq. ^ in units of nieC^, with pz = p± = as & 
function of z = z/£ for the Sauter potential V±{z) ((56|) for a = 5. 



A. WKB transmission probability for Klein-Gordon Field 



The probability of quantum tunneling in the z-direction is most easily studied for a scalar field 
which satisfies the Klein-Gordon equation 



ihdn + -AJz) 
c 



m^c > (/)(x) = 0, 



(5) 



where xq = ct. If there is only an electric field in the z-direction which varies only along z, 
we can choose a vector potential with the only nonzero component ([T|), and make the ansatz 
0(x) = e~ ' e^^^^-'-' 4'p^^£{z), with a fixed momentum p_L in the x,y direction and an energy £, 
and Eq. ^ becomes simply 

1 



-n'^+pl + my-^,[£-Viz)f 



By expressing the wave function (j)p^^e{z) as an exponential 



'Px,£ 



z) = 0. 



(6) 



(7) 



where C is some normalization constant, the wave equation becomes a Riccatti equation for S'p .£•: 



- ind'',Sp^,eiz) + [d,Sp^,e{z)f-pliz) = 0. 



(8) 



where the function Pz{z) is the solution of the equation 



pI(z) = \[£ -V{z)f -pl-my. 
The sohition of Eq. ([8]) can be found iteratively as an expansion in powers of fi: 
S^„e{^) = Sf^^,iz) - ihS^^leiz) + {-ihfS^^l,{z) + . . . . 



(9) 



(10) 
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Neglecting the expansion terms after Sl^^{z) = — log pz (z) leads to the WKB approximation 
for the wave functions of positive and negative energies can (see e.g. 30l. l31|) 
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^ J^n! £(^)/?i 
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P±.£^ 



(11) 



where 51 f(z) is the eikonal 



S'^lei^) = I pM)dz'. 



(12) 



Between the turning points z- < z < 2+, whose positions are illustrated in Fig. [H the momentum 
Pz{z) is imaginary and is useful to define the positive function 



1 



Kz{z) = dpi + mlc^ -^[£- V{z)Y > 0. 



(13) 



The tunneling wave function in this regime is the linear combination 
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2(«:.)V2 



exp 



K,dz 
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C 



■ exp 
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K^dz 



(14) 
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Outside the turning points, i.e., for z < 2;_ and z > z^, there exist negative-energy and positive- 
energy solutions for £ < £^ and £ > £-^- for positive pz- On the left-hand side of z_, the general 
solution is a linear combination of an incoming wave running to the right and outgoing wave 
running to the left: 



iPzY/' 



exp 



/M.]+^exp[-l/%.d. 



(15) 



On the right hand of z+, there is only an outgoing wave 



r 



(p,)l/2 

The connection equations can be solved by 



exp 



Pzdz 



2+ 



(16) 



C = 0, C± = e^'^/^C/2, T = C^ 



exp 



1 /■^+ 



K^dz 



(17) 



The incident flux density is 

J, = -^ [cP*dA - [d^^lc^] = ^^*<P = i^, (18) 

Zniei rrie me 

which can be written as 

Uz)=v,{z)n-iz), (19) 

where Vz{z) = pz{z)/me is the velocity and n^(z) = <j)*{z)(j){z) the density of the incoming particles. 
Note that the z-dependence of Vziz) and n^{z) cancel each other. By analogy, the outgoing flux 
density is \T\'^ /me- 

B. Rate of pair production 

Prom the above considerations we obtain for the transmission probability 

transmitted flux , , 

WwKB = . ., , n 20 

incident nux 
the simple exponential 

r 2 f^+ 1 

Wwkb{p±,£) = exp -- / Kzdz . (21) 

In order to derive from (j20p the total rate of pair production in the electric field we must multiply 
it with the incident particle flux density at the entrance z^ of the tunnel. The particle velocity at 
that point is Vz = d£/dpz, where the relation between £ and z_ is given by Eq. (j3|): 

This must be multiplied with the particle density which is given by the phase space density 
d?p/{2T:h)^. The incident flux density at the tunnel entrance is therefore 

3z U- )-^sj g .2vr/i)2 27rh ~ 'J 2^% (2^%)^ ' ^'^> 



m 



dpz (2'Khf 2-Kh J 27rh (27r/i)2 ' 

and the extra factor Ds is equal to 2 for electrons with two spin orientations 
It is useful to change the variable of integration from z to C(z) defined by 

and to introduce the notation for the electric field E{p±,£; C) = E[z{p±,£; (")], where z{p±,£; () is 
the inverse function of (j24p . the equations in ^ reduce to 

C-ip±,£;z.) = -l, C+ip^,£;z+) = +l. (25) 



In terms of the variable Ci the WKB transmission probabihty (j2ip can be rewritten as 



I 2?77. C 



1 + 



,2^4 



["'mSjE-l '^^) 



Here we have introduced a standard field strength £'0 to make the integral in the exponent dimen- 
sionless, which we abbreviate by 



Gip^,£)^- dC—^^^^—r. (27) 

vrj^i E{p_i_,£;C)/Eo 

The first term in the exponent of ()26p is equal to 2Ec/Eq, where 

Ee = mlc^/eh (28) 

is the critical field strength which creates a pair over two Compton wavelengths 2Ac = 2h/meC. 

At the semiclassical level, tunneling takes place only if the potential height is larger than 2rn-eC^ 
and for energies £ for which there are two real turning points z±. The total tunneling rate is 
obtained by integrating over all incoming momenta and the total area V± = J dxdy of the incoming 
flux. The WKB-rate per area is 

rwKB dJ ^ I 4%,W^Mp^,£). (29) 



Vi_ J 2TrhJ (27r^) 

Using the relation following from ([22 



d£ = eE{z_)dz_, (30) 

we obtain the alternative expression 



where <?(^;-) is obtained by solving the differential equation (j30p . 

The integral over p_|_ cannot be done exactly. At the semiclassical level, this is fortunately 
not necessary. Since E(. is proportional to l//i, the exponential in (j26p restricts the transverse 
momentum p_|_ to be small of the order of v^, so that the integral in (|3ip may be calculated from 
an expansion of G{p±,£) up to the order p\: 

2/''^^ yr^C^ r ldE{0,£X)/dC 



G{p^,£)^^J^^dC-^ 



■E{0,£;0/Eo 
where 6 = 5{p±) = {cp±)'^/{m'^c'^), and 



^ 2 E{0,£,C) '-'^^ 



G{0,£) + Gsi0,£)5 + ...i32) 



._lo,0,.,.i£.C^^5^^. (33) 



We can now perform the integral over p_|_ in (j3ip approximately as follows: 

'^ P^ .^-AEa/Eo){l+5)[G{0,e)+Gs(0,£)5] _ "^e^ f"^ ^^ ^-7,{E^/ Eo)[G{ld,e)+&G{Q,e) 






eEQ _-7t{e,/Eo)G{o,£) 



4TT^hcG{0,£) 



where 



G{0,£) = GiO,£)+GsiO,£). (35) 

The electric fields E{p±,£; Q at the tunnel entrance z- in the prefactor of pip can be expanded 
similarly to first order in 6. If z^ denotes the solutions of (j22p at p± = 0, we see that for small 6: 

A._ = ._-._«^^-. (36) 

so that 

E{z-)^E{z'^)-m,c'^^-. (37) 

Here the extra term proportional to 6 can be neglected in the semiclassical limit since it gives a 
contribution to the prefactor of the order h. Thus we obtain the WKB-rate pip of pair production 
per unit area 

rwKB ^ f^/.rwKB(^) ^ f^^ e^EoEiz) (^,/^„)G(0,^(,)) ^ . 

where z is short for z^. At this point it is useful to return from the integral / dz^eE{z^) introduced 
in (j3ip to the original energy integral / d£ in (j29p , so that the final result is 

TwKB _ f ^^ dsTy^KBJz) _^^ eEp r d£ I ^-^(EjEo)G(o,e) ^ ^39^) 



V± J V± UTT^hcJ 2TTnG{0,£) 

where ^^-integration is over all crossing energy- levels. 

These formula can be approximately applied to the 3-dimensional case of electric fields E(x, y, z) 
and potentials V{x, y, z) at the points (a;, y, z) where the tunneling length a = Zj^ — Z- \s much 
smaller than the variation lengths (5x_l of electric potentials V{x, y, z) in the xy-plane, 

a V ox_L 

At these points {x,y,z), we can arrange the tunneling path dz and momentum pz{x,y,z) in the 
direction of electric field, corresponding perpendicular area d'^V± = dxdy for incident flux and 
perpendicular momentum p^. It is then approximately reduced to a one-dimensional problem in 



the region of size 0{a) around these points. The surfaces z_(x_,y_,<f) and 2;+ = (x+,y+,<f) 
assciated with the classical turning points are determined by Eqs. ()24p and Eqs. ()25p for a given 
energy £. The WKB-rate of pair production (|38]l can then be expressed as an volume integral over 
the rate density per volume element 

TwKB = / dxdydz J^™ = / dtdxdydz ^^^ . (41) 

J ax ay dz J at ax ay az 

On the right-hand side we have found it useful to rewrite the rate Fwkb as the time derivative of 
the number of pair creation events dN^j^^iB/dt^ so that we obtain an event density in four-space 



C^^A^WKB ^2?, ^^^O^(^) ^-^{EJEa)G{0.e{z))^ (-42) 



dtdxdydz STT^nG{Q,£{z)) 

Here x,y and z are related by the function z = z-{x,y,£) which is obtained by solving (j30p . 

It is now useful to observe that the left-hand side of (|42p is a Lorentz-invariant quantity. In 
addition, it is symmetric under the exchange of time and z, and this symmetry will be exploited 
in the next section to relate pair production processes in a z-dependent electric field E{z) to those 
in a time-dependent field E{t). 

Attempts to go beyond the WKB results ([38]) or (p9]) require a great amount of work. Corrections 
will come from three sources: 

I from the higher terms of order in (/i)" with n > 1 in the the expansion (jlOp solving the 
Riccati equation ([8|). 



II from the perturbative evaluation of the integral over p± in Eqs. (j29p or (j3ip when going 
beyond the Gaussian approximation. 

Ill from perturbative corrections to the Gaussian energy integral (j39p or the corresponding 
2;-integral (i38]) . 



All these corrections contribute terms of higher order in h. 



C. Including a Smoothly Varying B(z)-Field Parallel to E(z) 

The above results can easily be extended to allow for the presence of a constant magnetic field 
B parallel to 'Ei(z). Then the wave function factorizes into a Landau state and a spinor function 

n 

first calculated by Sauter [1]. In the WKB approximation, the energy spectrum is still given by 
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Eq. ([2]), but the squared transverse momenta p\ is quantized and must be replaced by discrete 
values corresponding to the Landau energy levels. From the known nonrelativistic levels for the 
Hamiltonian p\/2me we extract immediately the replacements 

^2 



c^p^^ = 2meC^ X -^ — > 2meC^ x Tilol \n^ h^cr ) 



n = 0,1,2,- ••, (43) 



where g = 2 + a/vr + . . . is the anomalous magnetic moment of the electron, lol = eB/nieC the 
Landau frequency, with a = ±1/2 for spin-1/2 and cr = for spin-0, which are eigenvalues of the 
Pauli matrix cr^. The quantum number n characterizing the Landau levels counts the levels of 
the harmonic oscillations in the plane orthogonal to the z-direction. Apart from the replacement 
()43p . the WKB calculations remain the same. Thus we must only replace the integration over 
the transverse momenta J (Pp±/(2'irh)'^ in Eq. (j34p by the sum over all Landau levels with the 
degeneracy eB/{27rhc). Thus, the right-hand side becomes 

g-^ ^-7r(Ec/Eo)G(0,£) ST p-niB/Eo)in+l/2+ga)G{0,e) 
2TTnc 

The result is, for spin-0 and spin-1/2: 



^^-7,{B/Eo)in+l/2+ga)G{0,£)^ ^^^^ 



where 



2E0 .-AE./E,)GiO,e)f^^^^^(^s^^Q^^yE,) (45) 



foix) = —r , fi/2{x) = 2—- cosh — - (46) 

smhvrx ' smhvrx 2 



In the limit B ^ 0, Eq. dHj) reduces to Eq. (IMI) . 

The result remains approximately valid if the magnetic field has a smooth z-dependence varying 
little over a Compton wavelength Ag. 

In the following we shall focus only on nonuniform electric fields without a magnetic field. 

III. TIME-DEPENDENT ELECTRIC FIELDS 

The above semiclassical considerations can be applied with little change to the different physical 
situation in which the electric field along the z-direction depends only on time rather than z. Instead 
of the time t itself we shall prefer working with the zeroth length coordinate xq = ct, as usual in 
relativistic calculations. As an intermediate step consider for a moment a vector potential 

A^ = {Ao{z),0,0,A,{xo)), (47) 
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with the electric field 



E = -dzAo{z) - ao^2(xo), xo = ct. 



(48) 



The associated Klein-Gordon equation ([5]) reads 

i2 



ihdo + -Ao{z) 



+ ra 



2^2 



ihdz + -A^{xo) 



m^c > 4>{x) = 0. 



(49) 



The previous discussion was valid under the assumption Az{xq) = 0, in which case the ansatz 
4>{x) = e^*^*'^e*P^^^(/)p^^£-(2) led to the field equation ([6]). For the present discussion it is useful 
to write the ansatz as (/>(x) = e~^'P°^°'^e^'^^'^^'^(j)y,^^p^^{z) with pq = £ /c, and Eq. ^ in the form 



1 



£ — e dz E{z 



n2 ^2 ^ 



(50) 



Now we assume the electric field to depend only on xq = ct. Then the ansatz (t>{x) 
^tpzz/h^ip^:>ij_/h^^^^^^^^^-^ leads to the field equation 



-h' 



2„2 



d^-P^"^^' 



-Pz / dxQE{xQ 



n2 



'/'Px,p.(a^o) =0. 



(51) 



If we compare Eq. (|51|) with (|5U|) we realize that one arises from the other by interchanging 



z <^ xq, p± — > ip±, c ^ ic, £■ — > —iE. 



(52) 



With these exchanges we may easily calculate the decay rate of the vacuum caused by a time- 
dependent electric field -E(xo) using the above-derived formulas. 

IV. APPLICATIONS 



We now apply Formulas (j39p or (j38p to various external field configurations capable of producing 
electron-positron pairs. 



A. Step-like electric field 



First we check our result for the original case of a constant electric field E{z) = eEq where the 
potential energy is the linear function V{z) = —eE^z. Here the function (j27p becomes trivial 



G{0,£) = ^l\c^Jl-e = l, Gs{0,£) = 0, 



(53) 
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which is independent of £ (or z^). The WKB-rate for pair-production per unit time and volume 
is found from Eq. ([38]) to be 

1—1 

where V = dz^V\ . This expression contains the exponential ^-'^^cIEq found by Sauter [lil, and the 

Q Q 

correct prefactor as calculated by Heisenberg and Euler '2|, and by Schwinger [3]. 

In order to apply the translation table (j52p to obtain the analogous result for the constant 
electric field in time, we rewrite Eq. ()54p as 

where dN-\^KB/dxo = E^j^g/c and Ny/xB is the number of pairs produced. Applying the transla- 
tion table ([52|) to Eq. ([55]) . one obtains the same formula as Eq. ([M|) . 

B. Sauter electric field 

Let us now consider the nontrivial Sauter electric field localized within finite slab in the xy- 
plane with the width i in the z-direction. A field of this type can be produced, e.g., between two 
opposite charged conducting plates. The electric field E{z)z in the z-direction and the associated 
potential energy V{z) are given by 

E{z) = £;o/cosh2 ^^/^^ ^ y^^^ ^ -anieC^ tanh {z/£) , (56) 

where 

a = eEoi/meC^ = {£/Xc){Eo/E,). (57) 

Erom now on we shall use natural units in which energies are measured in units of m(.(? ■ Eigure[T] 
shows the positive and negative-energy spectra £±{z) of Eq. ([2]) for p^ = p_|_ = to show the 
energy-gap and energy-level crossings. Erom Eq. (JH) we find the classical turning points 



Tunneling is possible for all energies satisfying 

- \/rT5 + CT > £: > \fTVb - a, (59) 

for the strength parameter a > \/l + 5. 
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We may invert Eq. (j24p to find the relation between (^ and z: 



z = z{p± ,£;() = i arctanh 



s + CVT+s 



■In 



a+£+CVTTS 



a 2'" a-£-CVTT6' 

In terms of the function z{p±,£; ("), the equation ([^5|) . reads simply z± = z{p±,£] ±1) 
Inserting (j60p into the equation for E{z) in Eq. (j56p . we obtain 

2' 



(60) 



E{z) = Eo 



1 



CVTT6-£ 



a 



E{p±,£;C)- 



(61) 



We now calculate G{0,£) and Gs{0,£) of Eqs. 1^, ([MD and p3]l : 

o r o 1 1/2 r o 1 1/2 

G(0,f) = 2a^-a ia-£f-l - a {a + £f - 1 



(62) 



and 



O" 



GiO,£)+Gs{0,£) = :^\ ia-£f-l 



-1/2 



(a + f)2-l 



-1/2 



(63) 



Substituting the functions G{0,£) and Gs{0,£) into Eqs. (f38l) and (i39]) . we obtain the general 
expression for the pair-production rate per volume slice at a given tunnel entrance point z^{£) or 
the associated energy £{z^). The pair-production rate per area is obtained by integrating over all 
slices permitted by the energy inequality ([59]) . 

In Fig. [2] we show the slice dependence of the integrand in the tunneling rate (j38p for the Sauter 
potential ()56p and compare it with the constant-field expression ()54p of Euler and Heisenberg, if 
this is evaluated at the ^-dependent electric field E{z). This is done once as a function of the tunnel 
entrance point z and once as a function of the associated energy £. On each plot, the difference 
between the two curves illustrates the nonlocality of the tunneling process 33| . 

The integral is dominated by the region around £ ^ 0, where the tunneling length is shortest [see 
Fig. [1] and tunneling probability is largest. Both functions G{0,£) and Gs{0,£) have a symmetric 
peak at i5 = 0. Around the peak they can be expanded in powers of £ as 



G(0,£:)=2[(j2-(T(a2- 1)1/2] + 
and 
G{0,£)+Gs{0,£) 



(^2_''i)3/2 ^' + ^(^') = G-oC^) + lG2{a) £' + Oi£'), (64) 



(T 



1 (1 + 2C72 



+ - /. • ../,^ £^ + 0{£^) = Go{a) + -G2ia) £^ + 0{£^). (65) 



(a2- 1)1/2 2(c72-l)5/2 
The exponential Q^'^C!{o,t)Ec/Eo j^g^g g^ Gaussian peak around £ = whose width is of the order of 
1/Ec oc h. This implies that in the semiclassical limit, we may perform only a Gaussian integral 
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FIG. 2: We plot the slice dependence of the integrand in the tunneling rate ([38]) for the Sauter potential 
(1561) : left, as a function of the tunnel entrance z (compare with numeric results plotted in Fig. 1 of Ref. 
[26j): right, as a function of the associated energy £", which is normalized by the Euler-Heisenberg rate ((54|) . 
The dashed curve in left figure shows the Euler-Heisenberg expression (j54p evaluated for the z-dependent 
field E{z) to illustrate the nonlocality of the production rate. The dashed curve in right figure shows the 
Euler-Heisenberg expression (|54p is independent of energy-level crossing £. The dimensionless parameters 
are cr = 5, Eq/Ec — 1. 



and neglect the if^-dependence of the prefactor in (|39p. Recalling that £ in this section is in natural 
units with mf^c^ = 1, we must replace J d£ by nieC^ J d£ and can perform the integral over £ 
approximately as follows 



TwKB 

V± 



'D. 



eEonieC^ 1 



-^e 



-k{Ec/Eo)Go 



^TT^hc Gn 



d£ 
2^ 



^-n{E,/Eo)G'^ey2^j^ 



eEn 



-7t{Ec/Eo)Go 



(66) 



'iTr'hcGo2TTn^/G^Ej2Eo 

For convenience, we have extended the limits of integration over E from the interval (— 1 -|- a, 1 — a 
to (— oo,cxd). This introduces exponentially small errors and can be ignored. 
Using the relation (f57|) we may replace eEonieCp /he by e^E'^Hlja, and obtain 



TwKB [total] 
This approximate result agrees 



D 



e^El [^ ia'-ir/\ .^c,,^,EJE. 



"STT^h^cy Ec a5/2 



(67) 



34( 1 with that obtained before with a different, somewhat more 



23( 1 after the fluctuation determinant was 



complicated technique proposed by Dunne and Schubert 
calculated exactly in [24] with the help of the Gelfand-Yaglom method following Ref. [35|. The 
advantage of knowing the exact fluctuation determinant could not, however, be fully exploited 
since the remaining integral was calculated only in the saddle point approximation. The rate (j67p 



agrees with the leading term of the expansion (42) of Kim and Page |2ll |. Note that the higher 
expansion terms calculated by the latter authors do not yet lead to proper higher-order results 
since they are only of type II and III in the list after Eq. (j42p . The terms of equal order in h in 
the expansion (llOp of the solution of the Riccatti equation are still missing. 

Using the translation table (|52p . it is straightforward to obtain the pair-production rate of the 



Sauter-type of electric field depending on time rather than space. According to the translation 
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table (j52p . we have to replace i — > cT ^ where T is the characteristic time over which the electric 
field acts — the analog of i in (I56p . Thus the field ()56p becomes 

E{t) = ^o/cosh^ (t/T) , V{t) = -aruec^ tanh (t/T) . (68) 

According to the same table we must also replace a ^ ia, where 

a = eEoT /nieC. (69) 

This brings Gq{(t) of Eq. (j62p to the form 

G,{a) -. Gl{d) = 2[a(a2 _ i)V2 _ ^2]^ (70) 

and yields the pair-production rate 

rfvKB [total] ^ e^^o^ lEola' + l 



' '-" o,#!|i Jt f£!±iV^--i.<*)-./-o, (Ti) 



where F^yj^g [total] = 5Nwkb/5z is the number of pairs produced per unit thickness in a spatial 
shell parallel to the xy-plane. This agrees with Ref. J24i ]. 

Note also that the constant-field result (|54p of Euler and Heisenberg cannot be deduced from 
()67p by simply taking the limit £ ^ oo as one might have expected. The reason is that the saddle 
point approximation (I66p to the integral (I39p becomes invalid in this limit. Indeed, if £ oc o" is large 
in Eqs. ([62]) and ([63|) . these become 

G(0, £) ^ G(0, g) + G^(0, -g) ^ ^ _ L^2 ' (72) 

and the integral in (j39p becomes approximately 

^-niE,/Eo) r ^fi_ £2/a^) ^-HEc/Eo){£Va^) (73) 

For not too large i (x a, the integral can be evaluated in the leading Gaussian approximation 

_„o27r/i 27r/iV^/' ^ ^ 



corresponding to the previous result ()67p for large-cr. For a constant field, however, where the 
integrands becomes flat, the Gaussian approximation is no longer applicable. Instead we must first 



set o" — > oo in the integrand of (I73p . making it constant. Then the integral (I73p becomes J36] 



e-AEc/Eo)2fj/2T^n = e~'''-^^/'^°hieEo/meC^ 21x1%. (75) 



16 



Inserting this into (j39p we recover the constant-field result ()54p . We must replace 2£ by L to comply 
with the relation (I30p from which we obtain 

I d£ = f dzeE{z) = sEq f dz / cosh'^ {z / i) = 2ieEo = LcEq. 

In order to see the boundary effect on the pair-production rate, we close this section with a 
comparison between pair-production rates in the constant field (j54p and Sauter field (j67p for the 
same field strength Eq in the volume V±L The ratio -Rrate of pair-production rates (|67p and 
in the volume V\ I is defined as 



R 



rate 



'^O^nEJEg (O"^ - '^f'^ -wGnMEJEn 



Er. 



a 



5/2 



(76) 



The soft boundary of the Sauter field (j56p reduces its pair-production rate with respect to the 
pair-production rate ()54p computed in a constant field of width L = 2i. The reduction is shown 
quantitatively in Fig. [3l where curves are plotted for the rates (I54p and (j67p . and and for their 
ratio (|76p at £"0 = Ec and a = i/Xc [recall (j57p ] . We see that the reduction is significant if the 
width of the field slab shrinks to the size of a Compton wavelength Ac*. 
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FIG. 3: Left: Ratio i?rato defined in Eq. (|76p is plotted as function of a in the left figure. Right: Number of 
pairs created in slab of Compton width per area and time as functions of a. Upper curve is for the constant 
field ([5^ . lower for the Sauter field ([57]) ). Both plots are for Eq = Ec and a = £/Xc- 



C. Constant electric field for z > 

As a second application consider an electric field which is zero for z < and goes to —Eq over 
a distance i as follows: 



Eiz) 



Eq 
2 



tanh 



+ 1 



Viz) 



-—nipC < In cosh 



l)^ (.1' 



(77) 
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FIG. 4: Energies dH) for a soft electric field step E{z) of Eq. ([77]) and the potentials V±{z) ([771) for cr = 5. 
Positive and negative-energies £±{z) of Eq. ([2]) are plotted for p^ ^ p± = as functions of z = z/£. 

where a = eEoi/nieC^- In Fig. HI we have plotted the positive and negative-energy spectra £±{z) 
defined by Eq. ([2]) for pz = p± = to show energy gap and level crossing. From Eq. ^ we find 
now the classical turning points [instead of ([58]) ] 



Z4- = -In 
2 



"2g(£^±v^T+7)/a _ ^ 
For tunneling to take place, the energy £ has to satisfy 



£ < Vl + 5-0- In 2, 
and a must be larger than \/l + 5^. Expressing z/£ in terms of C as 



i 
z = z{pj_,£;C) = 2^" 



2g{^+Cv'i+3)/<7 _ ^ 



so that z± = z{p±,£; ±1), we find the electric field in the form 



E{z) = Eo 



1 



1 



{(VT+S-£)/^ 



^Eip^,£;C)- 



Inserting this into Eq. ([27|) and expanding Eq/E{p±,£-X) in powers we obtain 

g-nf/cr 2 rl 

n=l 



(78) 



(79) 



(80) 



(81) 



G{p^,£) = l + J2 



^ 2" vr 7-1 

n=l 



dCVl - C' e"^/" = 1 + E e-"^/"/i(nyrT5/a), (82) 



where /i(x) is a modified Bessel function. Expanding Ii{n\/1 + (^/o") in powers of 5: 



/i(nyrT5/a) = Ii{n/a) + (n/4a)[/o(n/a) + /2(n/a)](5 + 



(83) 
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we identify 

oo 

G(0,£:) = l + 5]e-"^/'^/i(n/a), (84) 

n=l 
-. oo 

G{0,£) + Gs{0,£) = l + -^e-"^/'^[(n/a)Io(n/a)-/i(n/a)]. (85) 

The integral over £ in Eq. (j39p starts at iS< = 1 — o" log 2 where the integrand rises from to 1 as 
£ exceeds a few units of a. The derivative of e^'^^^' '''^^' ' drops from 1 to e^'^^ "' '^' over this 
interval. Hence the derivative d£e~'^''^'''^°''^^^'^' is peaked around some value £. Thus we perform 
the integral / d£e~'^^^''/^°'^^''^'^^ by parts as 

The first term can be rewritten with the help of d£ = eE^dz as e~'"^ "' '-'''\eEQ\l /2, thus giving rise 
to the decay rate (|54p in the volume Vj_^/2 , and the second term gives only a small correction to 
this. The second term in Eq. (j86p shows that the boundary effects reduce the pair-production rate 
compared with the pair-production rate (|54p in the constant field without any boundary. 

V. TUNNELING INTO BOUND STATES 

We turn now to the case in which instead of an outgoing wave as given ()16|) there is a bound 
state. We consider a linearly rising electric field whose potential is harmonic: 

.(.,.. „(^), .,.)^^(^)^ (S7) 

It will be convenient to parametrize the field strength Eq in terms of a dimensionless reduced 
electric field e as -Eq = e/ic/eA^ = eEc- In Fig. [S]we have plotted the positive and negative-energy 
spectra £±{z) defined by Eq. ([2]) for p^ = p_|_ = to show energy gap and level crossing for e > 
(left) and e < (right). If e is positive, Eq. dH) yields for z > 0, 

z± = \c\j-(£^VlTtf''^ , z+<z^, (88) 

and mirror-refiected turning points for z < 0, obtained by exchanging z± -^ —z± in (j88p . Negative- 
energy electrons tunnel into the potential well —z^ < z < -|-2;+, where <?><?+, forming bound 
states. The associated positrons run off^ to infinity. 
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FIG. 5: Positive- and negative-energy spectra £±{z) of Eq. ^ for p^ = p^ = as a function of z = z/Xc for 
the linearly rising electric field E{z) with the harmonic potential ((87)) . The reduced field strengths are e — 2 
(left figure) and e — —2 (right figure). On the left, bound states are filled and positrons escape to z = ±oo. 
On the right, bound electrons with negative energy tunnel out of the well and escape with increasing energy 

to 2 = ±00. 

A. WKB transmission probability 



Due to the physical application to be discussed in the next section, we shall study here only 
the tunneling process for e > on the left-hand side of Fig. [5l We consider the regime z < 
with the turning pints — Z- < — z^. The incident wave and flux for z < — z_ pointing in the 
positive z-direction are given by Eqs. (J15p and (jlSp . The wave function for —z^ < z < —z^ has 
the form Eq. ()14p with the replacement z^ -^ — z_. The transmitted wave is now no longer freely 
propagating as in p6|) . but a describes a bound state of a positive-energy electron: 



^sjz) 



B 



■ cos 



Pzdz 



z+ 



(89) 



(p,)l/2 

The Sommerfeld quantization condition 

1 /•+^+ 

fixes the energies £n- The connection rules for the wave functions (J14p and ()89p at the turning 
point —Zj^ determine 



-/ p^dz = 7r(n-F i), n = 0,1,2,, 



(90) 



B = V2C+e- 



exp 



K,dz 



(91) 
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Assuming the states (j)£^ {z) to be initially unoccupied, the transmitted flux to these states at the 
classical turning point — 2_|_ is 



h 



nil 



<t)sSz)dAlM 



\B\ 



\C, 



exp 



-z+ 



Hi,dz 



(92) 



z^-z+ (2me) rrif. 

From Eqs. (jlSp . (j92p . and (j20p we then find the WKB transmission probability for positrons to fill 
these bound states leaving a positron outside: 

2 r-^+ 



W^KB{p±,£n) = exp 



n 



K^dz 



(93) 



which has the same form as Eq. ()21|] . The same result is obtained once more for z > with turning 
points 2+ < z_, which can be obtained from (j93p by the mirror reflection —z± <-> z±. 



B. Energy spectrum of bound states 

From Eq. Q for pz and Eq. ([87]) for the potential V{z), we calculate the eikonal ([90]) to 
determine the energy spectrum £n of bound states 



1 /•+^+ 
h J-Z+ 



Pzdz 



'A3 



C7^0 



^+ 



(z2-4)(z2-£) 



1/2 



dz 



2€Z+ 

3A^ 



(4 + £)i5;(t)-(z^-4)ir(t) 



(94) 



where E{t), K{t) are complete elliptical integrals of the first and second kind, respectively, and 
t = 2;+/z_. The Sommerfeld quantization rule (f90ll becomes 



2{en-Vi + 5) 



1/2 



8nE{tn) - [Vl + ~5)K{tn, 



1/2 



"^" + 2^' ^"= Un + 7rn -^''^ 



For any given transverse momentum p_\_ = v5, this determines the discrete energies £n- 



C. Rate of pair production 



By analogy with Eqs. ()29p and (I39p . the transmission probability (I93p must now be integrated 
over all incident particles with the flux (|23p to yield the rate of pair production: 

TwKB „ n sr^^n f 'f'Pl. 



V^ 



^"•?S/ 



(2,rR)2''''™"<=<^^'^"'' 



9 n 1*^^01 Y^ ^-^ 

^^' ■ 2ti^Z^ O-^' 



1 



,-7r(£,/£;o)G{o,£:„) 



(96) 
(97) 



■47r27icV2vrG(0,£:„) + G5(0,f„) 
In obtaining these expressions we have used the energy conservation law to perform the integral 
over £. This receives only contributions for S = £n where J d£ = ujrfi = £n — £n-i- The factor 2 
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accounts for the equal contributions from the two regimes z > and z < 0. The previous relation 
([301) is now replaced by 



1/2 



ujnh = \eE{z1)\Az'l. 
Using Eq. (j24p and expressing z/Xc > in terms of (^ as 

Z = z{p±,£n;C) = \c\l- \£n - CVl + dj 

we calculate z± = z{p±,£n', if)) ^-nd find the electric field in the form 

E{z) = Eo^^ [Sn - CVT+6y^^ = E{p^,£n; ()■ 

Inserting this into Eq. (j27p we obtain 

2 /? rl 



(98) 



(99) 



(100) 



G{p±,£n 



dC 



VT^C' 



7rV2 7^i ^[f„_cVrT5]i/2' 



8 f7{£i + l)y^ 



SvrV 2 (1 + <5)V4 



Srj^rJ-^ 



[l-£i)K{qt,)+£'^E{ql) 



(101) 



where £^^ = <S„/(1 + (5)"*^'^ and g„ = J2/{£f^ + 1). Expanding G{p±,£n) in powers of (5 we find the 
zeroth order term 



G{0,£n) = —J-{£n + l)'^^[{l-£n)K{qn)+£nE{qn)] 

ovr V ^ 



(102) 



and the derivative 

Veq-n 



Gs{0,£„ 



Svr £„ 



qn-l {^-hqn+£n{T-Qqn))E{qn) + {l-£n-^£l){qn-l)K{qn) .(103) 



where qn = \J1l{£n + 1)- 



VI. COULOMB ELECTRIC FIELD 



We now come to the physically interesting system of a bare nucleus of high charge Z in the 
vacuum which fills its empty bound states by quantum tunneling of electrons from the negative- 
energy continuum around it. The radial Coulomb field and the potential are given by 



^, , dihc ,,, , dihc 
eE{r) = ^, U(r) = , 



where a = Za. The Schrodinger equation of this problem reads 

2 



s + ^] +c^n'x/^-my 



#(x) = 0. 



(104) 



(105) 
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After factorizing out spherical harmonics V'£'(x) = -R£-^/(r)yjm(x), the radial wave functions satisfy 



s^ + ch^ + c^n^ 



c)t 



i{i + i) 



a 



m^c ) rR£^i{r) = 0, 



(106) 



where / = 0, 1,2,3, . . . are the quantum numbers of angular momentum. The term in brackets 
cannot be treated semiclassically since it contains a factor fi . This destroys the possibility of a 
systematic expansion of the wave function in powers of h. In fact, the % /r^-potential should be 
considered as being part of the kinetic term, not of the potential. In fact, both have the same 



scaling dimension. There exists a heuristic way of accounting for this due to Langer 



m- He 



introduced a change of coordinates to .^ = logr, so that 9^ rRg^i{i 



-2? 



{dl-d^)r{ORe,MO) 



-2? 



(5| - i)e-«/2r 



)e ^''^r{^)R£^i(r{£^)), which brings (11061) to the form 



rHC){£' - my) + ^chSa r(0 + c^h' {a| - [l{l + 1) + i - a^]} J ^)ReMO) = 0, (107) 
Now the h /r^-part of the kinetic term has become a trivial constant which no longer influences 



the semiclassical treatment 



m- 



Equivalently, we can apply a corrected semiclassical treatment directly to Eq. (jl06p if we add j to 
/(/ + 1) in the numerator of the centrifugal barrier. This so-called Langer correction is implemented 
in ()106|) by replacing l{l + I) — 6? by A;(Ai + 1), where 



A/ 



G+i; 



a 



1/2 



/ 



a 



2^ + 1 



+ 0(a^ 



(108) 



Note that Eq. (jl07p has lost the singularity of Eq. (jl06p at the origin, a fact which is often 
considered to be the motivation for going to the variable ^. However, due to the equal scale of 
fi jr'^ and the gradient term ti 9^, this is not the relevant property. The fi /4r ^-correction is 
needed at any distance to obtain the correct wave function at the semiclassical level. In fact, this 
wave function, although being approximate, turns out to produces the exact energy levels of the 
relativistic Coulomb system. 



1. Semiclassical quantization for point-like nucleus 



The semiclassical treatment of Eq. (|106p with / replaced by A^ starts out with the Langer- 
corrected energy 



£±iPr, I; r) = ±^{cVrY + {hcY{l + i)2 A^ + mlc^ + y{r). 
We now impose the Sommerfeld quantization rule upon the eikonal: 

S{£) = drpr = nh{nr + 5), 



(109) 



(110) 
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where n,. = 0, 1, 2, ... is the radial quantum number, and 

1 



Pr 



f 2 + 2chaS/r - c^n^{\i + if/r'^ - mld^, 



whose zeros yield the turning points for each /: 

ch 



a8 ± Va2£:2 + ^Xi + |)2(£:2 _ mlc'^ 



Consider first the energy regime < £ < nieC^- For < n < Tq, we rewrite pr as 

Pr 



Or = -^Jmld^ — 8'^ -^ {tq — r) (r — n), 



and perform the integral in (jllU|) to find the eikonal 



f 



71" , 



S{£) = -Jmlc^ - £^ -(ro + n - 2^/rvT) = vr/i 



a£ 



^mjc^ - f 2 



(A^ 



where 



i\2 



r^r; 



c27i2i4i±i)l, r, + n = cn- 



2a£ 



(111) 



(112) 



(113) 



(114) 



(115) 



m'^c'^ — i?2 ' m-gC^ — £"2 

Inserting S{£) into the quantization condition (jllOp . we obtain the following equation for the exact 

bound state energies 

1 



2mld^ 



o? £nl' 



2 mld^ {rir + \i + 1)2 ' 



(116) 



This is solved by 



■'nl 



± 



rrif-c 



^l + {ZaY/{n-5iY 
6? 



±mpC 



1 



--2 3a^ 



Q 



+ 0{a 



^fi^ 



2n2 8 n^ n^{2l + 1) 
where n = n^ + Z + 1 is the principal quantum number of the atom, and 



6i = (/ + i)-(Az + i) 



a 



21 + 1 



+ 0{a^). 



(117) 



(118) 



Although derived by a semiclassical approximation, these happen to be the exact values, due to 
the Langer correction in (jllip . 

For a Dirac particle, the energy-spectrum £nj can be obtained from a slight modification of 
Eq. (fTTTl) : 



S, 



nipC 



nj 



(119) 



^l + {Zay/in-5,r 

where j = 5, |, • • •,n — I is the total angular momentum. For each j, there are two degenerate 
states with orbital angular momentum I = l± = j it | , for which we define Aj^ =l± — Sj , with 

Sj^j + h- JiJ + \Y - «' = «V(2j + 1) + 0{a% (120) 



The radial quantum number n,. is related to the others by rir + A,_,_ + 1 = n^ + /± + 1 — 5; = n — 5j. 
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2. Semiclassical quantization for finite-size nucleus 



According to Eq. ()120p . the energy ()117p becomes imaginary for the I = 0-state of spin-0 particle 
when a > 1/2, and of a spin-| particle, when a > 1. This signalizes the crash of the electron into 
the Coulomb potential of a point-like positive nuclear charge. The result is, however, unphysical 
since nuclei always have a finite radius r^. We account for this by the approximation a uniform 



distribution of charge inside the nucleus 



V{r) 



7|, so that the potential is 
hca/r, r > r^, 



fl21) 



where /(x) = (3 — x^)/2. Now the Sommerfeld quantization condition for the eikonal (jllOp reads 



S{£) = S^^\£) + 5(2)(f ) = r drpi^^ + r drpf^ = 7rh{nr + i), 

Jfo Jrn 

with the radial momenta in the different regions: 



(122) 



P 



(2) 



f + ^c^/ - 
r„ \r 



n \' n 



c^h^^^ + ^ 



^^ -mlc^, 



£ + hc 



a 



n2 



cH^(l±M- 



rald^. 



(123) 
(124) 



Let us first calculate the eikonal S^^'{£) inside the nucleus. Expanding p)- in powers of r/r^ 
keeping terms up to the order 0[(r/rn)^] we obtain 



P^ 



^ 3(a 
S + cfi^ 
c \\ 2r 



2 4 " /c , „t.3QWr\2 (/+i)2 



mzc \£ + ch-— 

r„ V 2r 



11/ \ ' n 



1/2 



c'h' 



3a 



^ + c^§r) -m^ 



1/2 



^/ 2 ~2\l/2 
-(r -To) / 



(125) 



where 



fo = c?i 



(l + h)' 



{£ + ch§-J-mlc^ 



1/2 



3a 



(/ + i)r„. 



(126) 



The approximation is good for r^/Xc <C 1, which is assured as long as 



\£\ <C tic—, and Af = -(A/ + i)^ = a^ _ (/ + i)2 ^ q_ 



(127) 



The second inequality determines a maximum value of the angular- momentum I for a given a. 
Inserting the approximate value (jl26p into (jl25p , and the associated pr into (jl22p , we obtain the 
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r < Tn-part of the eikonal 



sW(f) 



TTlpC 



(r — Vq) ' + rn(/ + 5) arcsin 



ro 



^'^ 



21 + 1 
3a 



1/2 



+ {1 + 1 



ro 



arcsm 



2/ + l\ vr 
3a 



(128) 



We now calculate the eikonal S^'^\£) outside the nucleus in Eq. (jl22p . We begin with the 
negative-energy regime —m^c? < iS < 0, where Eq. (11270 is satisfied so that Eq. ()112p gives 
fo ^ i^n, and unphysical turning points n < and |ri| > Vq. Using Eq. (jllSp for pr '(r) and 
integrating it over r^ < r < r^ we find 

-roTi) ' In ■ 



R 



2{-ror\) + {ro+n)r + 2{-ronY/'^R ro+ n . {vo+n) -2r 



■ arcsm ■ 



(ro-n) 



(129) 



where R{r) = \/{ro — r) {r — n), and 



To — Ti = 2c^ 



a2f2 + (^^+i)2(^2_j„2c4) 



rn-Jc 



2.4 _ £:2 



Under the conditions r^/Xc <^ 1 and rn/|ro,j| <C 1, Eq. (|129p becomes approximately 

a£ 



s^^H£)^n 



, — __; ^ , h arcsm — — AH 1 + m — -pr 

^m'ic'^-£-^ \2 XI J V c^^z 



(130) 



(131) 



where 



XI = V^?("ieC^ - ^^) + a^£^- (132) 

For Tji <C Ac, the first eikonal part S^^'{£) of Eq. (jl28p is negligible compared with the second 
part S^'^'{£) of Eq. (I13ip . Thus can simply do the calculation with a Coulomb potential cut off at 
r = Tn- The charge distribution inside nucleus is irrelevant here. 
For negative energies £ close to zero, we approximate 
1 



s^^\£)^n 



2A; V V V?7ip.c2 



£ V Tia r £ 
+ 



2 ynipC^ 
and the quantization rule ()122p yields the energy levels 
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-vripC 



iraAi + W(7raA/)2 + 8(a2 + ^2^, ^ ^^^ + i) + A,(l + In 



AcAi- 



2(a2 + A2) 
For negative energies £ close to —meC^, expression (J13ip simplifies to 

'A^ Af 



S(2)(f)«7j 



2 6a2i^^-^'P + l-A^ 



(134) 



(135) 
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where rf = (m^c^ — iS^)/(?n,eC^)^, so that the quantization rule (jl22p gives 

7rK + i)+Az(2 + lnf?f,) 

£nl ~ -m-eC^ 



7rK + i)+Az 2 + ln^ 



Sa^Ai - Af 

The critical value cic(0 can be obtained from Eq. p22|) . and from (|135|) for 77^ = 0, 

2 



(136) 



2 ^ „2 



Af = di-{1 + 1)^ « vr 



2 + 1-]^ 



(137) 



Consider now the negative-energy regime below —rrif.cp' under the assumption (jl27p . so that 
Eq. (J112p yields Tq ^ rn, and unphysical zeros at < Tq < n. By writing pr in Eq. (jllip by 
analogy with in ()113p as 



Pr = -y -^^ - "iec"^ ~v ('^o - J^) (n - ?•), 



(138) 



we can perform the second integral in ()122p and obtain 
5(2)(£:)= r drpr = -j£^ 

Jrr, C 
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fy I ^\|2^ 2{rori)-{ro + ri)r + 2{roriy/^R r^ + r, = 

R — (ron) ' m m\2R + 2r — [ro + r-A 

r 2 



, (139) 



where R{r) = \/{ro — r) (r; — r). Under the conditions r^^/Xc <^ 1 and r^/ro^i ^ 1, we obtain 

rnXl 
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(140) 



Assuming hca/r^ ^ \£\^ rn^.c?, we expand S^'^'{£) in terms of me(? /\£\ and the leading order is 



and the quantization rule ()122p yields energy-levels: 



mec^ a /A; A rn 2 a2n 
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(141) 
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(142) 



As for a spin-| particle, the critical value ac[j) and energy-spectra £nj are obtained from the 
critical value (jl37p and energy-spectra (jl36p and ()142p by the following replacement: 



A2 = a2-(/ + i)2 ^ A^ = a2-(j + i)^ 



(143) 



As before, for each j there are two degenerate levels of orbit angular momentum I = l± = j it |. 
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3. WKB transmission probability 



The Schrodinger equation Eq. ()106p can be written as 



c'n^^ + cVAr) 



rRe,i{r) = 0, 



where 



^pl = [8-V[r)f-^n'^l±fi. 



2 4 



(144) 



(145) 



which looks hke the Klein-Gordon equation ([6]) in one dimension along the r-axis. As before, we 
assume the nuclear radius rn to be much smaller than the Compton wavelength <C Ac, Outside 
the nucleus, V{r) is the Coulomb potential. From th condition p^ = we calculate the classical 
turning points rn ^ r+ < r_, by analogy with Eq. ^. For a give energy 8 < —nieC^, we have 
three regions (see Fig. [6]) : 

(i) r- < r and £ < £- , 

(ii) r_|_ < r < r_ and £^ < £ < £^, 

(iii) r < r_|_ and £ > £jf-. 



Energy 




FIG. 6: For a radial Coulomb field, the positive and negative-energy spectra £±{r) of Eq. ^ are plotted as 
a function of f = r/Xc- They are found by solving condition pr — a,t I — and a — 1.27 {Z = 174). 



Starting (i) where p^ > 0, Eq. ()144p has two independent solutions corresponding to pr < 
(incident wave) and pr > (reflected wave) , 
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(146) 
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where C^ = e^*'^'^C''/2. The corresponding solution in the region (ii) is given by 
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(147) 



2(k^)V2 r 

where k^ = —ipr and p^ < 0. For a purely incident wave, only the second term in R^- i is present, 
with an incident flux density at r = r_: 
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(148) 



where (pr = fRg i and dr = {l/r){d/dr)r. The superscript - indicates the inward flux. 

We consider bound states of energy £ < —m^^c^ confined within the region (fo, r_|.), where fo < r^ 
is another classical turning point of positive-energy branch £j^ inside the nucleus, in addition to 
r_|_. Its value is obtained from Eq. (jl45p for p^ = 0: 



£ - V{fo) = +Jc^h^^^^J^ + m2c4 



^2 



(149) 



By analogy with Eq. (j89p . the semi-classical wave function of bound states is, 
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■f-^ = J7^-^o^ 



^r I 
Pr 

satisfying Sommerfeld's quantization rule 
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(150) 



n 



Prdr = vr {rir + 1) , n^ = 0, 1, 2, 



(151) 



For each n^ = 0, 1, 1 . . . , the solutions yield the discrete energies £ni with the principal quantum 
number n = rir + I + 1, where the angular momentum / can take the values I = 0, . . . ,n. The 
continuity of wave functions ()147p and (jl50p at the classical point r_|_ leads to 
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Assuming the bound states R^^^ be unoccupied, the transmitted flux to these states at the classical 
turning point r^ is 
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(153) 



From Eqs. (|148p . (jl53p . and (j20p we obtain the transmission probability of an electron to tunnel 
into the bound state: 

^2 " 
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(154) 
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Normalizing the incident flux density (jl48p at r^ we have 

(2/ + 1) /■ dpr 



I I 



Vr[r^, 



AtttI J 2-kV 
and the rate of pair-production in the state with angular momentum I becomes: 

rwKB('?,0 



v^ 



n ( a(2[+1) f dp, 
DsVr{r_)—;—-p^— I 7^::^exp 



4:7rri 



2Trh 



Krdr 



r+ 



D. 



{21 + 1) 

47rri 



211% 



exp 



7i 



Krdr 



r+ 



(155) 



(156) 



where t'r(r_) = d£ /dpr\r=r-- This is evaluated further in the same way as Eqs. (j96p and (j97p : 
the integral over £ has only contributions from the bound state energies £ = £nu so that j d£ is 
equal to (-ij„/?i, where lij„/ = £ni/'h is the frequency of the bound state. As a result, the sum over all 
Eq. ()156p takes the form 
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^.E^^Z^^exp 
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(157) 



-^. Sauter exponential factor in Coulomb potential 



For brevity we use natural units where r is measured in units of the Compton wavelengths 
^C = %/mi.c, and £ni in units of m^.c^. By setting pr = in Eq. (jllip . we obtain the analog of 
Eq. (H 



1 



£nia+cJ{i + m£ii-e)+^' 



^ 9 



(158) 



r{i,£nf,o {i + ^j^e 

The role of p\ is now played by (/ + i)^/r^. The classical turning point r^{£ni,l) where the 
negative-energy states can tunnel to the positive ones at r+(iS„;,/) is obtained from Eq. (jl58p by 
inserting (^ = ±1. We identify r+ = Tq in previous sections. Since r_ > r+, the electrons move 
inwards to neutralize the Coulomb potential, the positrons are pushed outwards to infinity. At the 
highest energy £ni = —1 where tunneling can occur, we find r_ = +cxd and r+ = [a^ — {1 + l)^]/^^- 
For each angular momentum / and energy-level £ni, we obtain the Sauter exponential factor 
from the analogs of Eqs. (f26l) and (i27|) : 
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and 
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(160) 



r^{l,£nr,C)_ 

Note that the bracket in the integral is the analog of the prefactor [l + {cp±)'^/m1c'^] in (j26p . Here 
it can no longer be taken out of the integral. Instead, there is now another simplification. Since 
E = ahc/r'^, the function G{l,£ni) becomes 

2 '■I 



G{l,£nl) 
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and with r{l,£;Q from ()158p : 
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(/+i)n.(i62) 



[(fn;a)2-C2((/ + l)2(£:2^-C2)+d2)]2 

The result of this integral is surprisingly simple remembering that r± are given by r±{£ni',l) of 
Eq. (HMD for C = ±1: 



Gil,£ni) = 2{l + \Y[£ii - \£niW£'^i - 1 +2d 




1 



(163) 



valid for £ni < —1 and a > dc(0- Inserting this into the Sauter exponential factor (jl59p . and this 
further into (|157p , we obtain the rate of filling the empty bound state levels around the Coulomb 
potential of the nucleus. The sum over I extends to the largest value permitted by Eq. ()127|) and 
dc of (fT37p . 

The result for a spin- 5 particle can be obtained by replacing £ni =^ £nj in Eqs. (|158p and (|163p . 
and 



(/ + i)2^(j + i)(j±l + l), 



(164) 



where j ± 1 respectively corresponds the state of parity (— )^ with orbital angular momentum 



VII. SUMMARY AND REMARKS 



By studying the process of electron-positron pair production from the vacuum by a nonuniform 
electric field as a quantum tunneling phenomenon we have derived in semiclassical approximation 
the general rate formulas ()38p and (I39p . They consist of a Boltzmann-like tunneling exponential, 
and a pre-exponential factor, and are applicable to any system where the field strength points 
mainly in one direction and varies only along this direction. The formulas require the evaluation of 
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the functions G{0,£) of Eq. ([32]) and Gs{0,£) of Eq. ([33]) . This has been done for several different 
field configurations. 

For electrons arriving by tunneling in a bound state of a harmonic electric potential, the general 
expressions are given by Eqs. ([97|) and (|157p as functions of the frequencies Un of the bound states. 
The discrete energy levels £n at fixed p± are found from the eikonal and Sommerfeld quantization. 
For fermions, the expression for pair-production rate should be multiplied by the Pauli-blo eking 
factor for the rate which makes it zero if the bound state £n is occupied. 

In the Coulomb electric field E(r) = eZ/r^ of a nucleus with finite radius rn <C Ac, we have 
given first the semiclassical energy- levels £ni in Eqs. ()117p . ()134p . (|136[) . ()142p . and the formulas 
for the associated pair production rate in (jl57p . (jl59p . (jl60p for £^1 < —rrieC^. The critical value 
etc = ZcO is found from Eq. (jl37p as a function of the principal quantum number n and the angular 
momentum I, which agree with the one found in Refs. [Tl, 139(1 the n = 1-state ISi . 

2 

The number of energy-levels £ni < —meC^ that are able to accommodate electrons produced 
from the vacuum is limited, and pair production ceases when all these levels are fully occupied 
even if the electric field is overcritical. 
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